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Abstract 

Calculation of the distribution of the average value of a Gaussian random field in a finite domain 
is carried out for different cases. The results of the calculation demonstrate a strong dependence 
of the width of the distribution on the spatial correlations of the field. Comparison with the 
simulation results for the distribution of the size of the cluster indicates that the distribution of 
an average field could serve as a useful tool for the estimation of the asymptotic behavior of the 
distribution of the size of the clusters for "deep" clusters where value of the field on each site is 
much greater than the rms disorder. 
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I. INTRODUCTION 



A common feature of any random medium is the formation of clusters. Probably, the 
most well-known problem where statistics of clusters has been extensively studied is the 
famous percolation problem l|. In the lattice percolation model any site is occupied with 
the probability p and non-occupied with the probability 1 — p, and each cluster is a set of 
connected occupied sites. An important characteristic of the random medium is distribution 
of the size of the clusters, or, more exactly, the number of clusters with s sites per lattice 
site, Ug. Knowledge of the statistics of random clusters is vital for the description of many 
important natural phenomena, such as the conductivity of disordered materials, the flow of 
liquids in porous media, fracture processes in materials, or even the dynamics of landscapes 
and forest fires . Exact or reliable approximate analytic results for Us are not 

very numerous. Possible examples include the number of percolation clusters for s ^ 1 and 
p^O 

Hs OC S~^p~\ (1) 

or the corresponding distribution near the percolation threshold Pc 



Us (X s exp [C{p - Pc)s^] , p Pc, s > 1. (2) 

Here 6, r, C, and 5 are some constants [3, 1^, Iq]. Most results in this area were obtained 
using scaling arguments with subsequent testing of their validity with extensive computer 



simulation 



Id, 



The reason for the scarcity of analytical results is obvious: it is 



difficult to take into account various shapes of clusters. In addition, most known results 
for Ug are obtained for the uncorrelated case, i.e. the case when lattice sites are occupied 
independently of each other. 

Cluster numbers for the case, where sites are occupied not independently, have been 

nnnnn 

studied for the problem of correlated percolation |12l . 1131 . 1141 . 1151 . Ilq . Correlation is usually 



introduced by the short range interaction between different sites (most popular cases are 



the Ising model 



17 



181 ] and the q-states Potts model 



Attention in this area was 



almost exclusively focused on the behavior in the vicinity of the percolation threshold. 

In this paper we are going to consider the statistics of Ug for another case of correlated 
random distributions, namely for a Gaussian random field U{r) (we will call U{r) the random 
energy for reasons that will be obvious later), where the cluster may be defined as a connected 
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set of sites, all of them having an energy greater than the threshold value Uq > 0; for a 
Gaussian random field this is equivalent to the cluster with sites, having the energy less 
than —Uq. We are going to consider the distribution of "deep" clusters with Uq ^ a (where 
a is the rms disorder and where we assume zero average for U), i.e. the situation far away 
from the percolation threshold. Most attention will be paid to the particular kind of a 
Gaussian random field having binary correlation function 

C{f) = {U{r)U{0)) ^Aa'-, r ^ a, a' = {U\f)) , (3) 



where the angular brackets denote a statistical averaging and a is the lattice scale. 



'his 



particular correlation function naturally arises in the model of dipolar glass (DG) [2l|, |22| . 
which is popular for the description of the charge transport properties of organic materials. 
In the simplest realization of the DG model we assume a random and independent orientation 
of dipoles occupying sites of a regular lattice, while charge carriers interact with dipoles by 
the long range charge-dipole interaction. In this model the energy of a charge carrier is 

— e r — r„ 

where d is the dipole moment of the molecule, and e is the dielectric constant. Using an exact 
analytic calculation as well as computer simulations it was shown that for the DG model 
random energy U (r) is a Gaussian random field if the average distance between dipoles is 



not significantly greater than the lattice scale 23|, |2J, |25|]. The correlation function of the 



DG has the form ([3]) and in the case of a simple cubic lattice A ^ 0.76 [26|]. This model 
was suggested to explain the Poole-Frenkel dependence of the carrier drift mobility fi in 
polar disordered organic materials on the applied electric field E over a broad range of field 
strengths 



In /i oc V -E (5) 
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29j. A power law decay of the correlation function ([3j) means an extremely 



long range correlation in the random energy landscape in organic materials. For this reason 
clusters have wide size distribution (see Fig. [T]). 

In this paper we show how to calculate analytically another characteristic of the correlated 
medium, namely, the distribution of the average random energy in a domain. An attractive 
feature of this distribution is that it is much more easy to calculate. We argue, then, that this 
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FIG. 1: Distribution of site energies U in the lattice model of dipolar glass. A sample with the size 
of 50 X 50 X 50 lattice sites is shown. Black and white spheres represent the sites with positive and 
negative values of C/, correspondingly, while the radius of a sphere is proportional to the absolute 
value of U. Sites with small absolute values of \U\ (less than a) are not shown for the sake of 
clarity. 

distribution provides valuable information about the cluster numbers for large clusters and 
Uo/a ^ 1, that is in the region far away from the percolation threshold. Our consideration 
will be limited to the 3D case, though generalization to other dimensions is obvious. 

II. DISTRIBUTION OF AN AVERAGE VALUE OF THE RANDOM FIELD IN A 
DOMAIN 

Let us calculate the distribution Pv{Uo) of the average value Uo of the random energy 
f/(r) in a domain with volume V (here we consider a spatial average and use the same 
notation Uq for the average energy). Let us start with a continuous model of the medium, 
which is valid for s ^ 1. In this model the distribution Py{Uo) is the average of the delta 
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function 

PviUo) = {^{yj dfV{r}Mf^ " ^0 ) ' 
where equals 1 inside the domain and outside, and may be presented as a path 

integral over all realizations of the scalar field U (r) 

PviUo) = ^j'^U6(^j dTU{T)fy{T) - e-^ (7) 

Z = j VUe-^, drdfiU{f)G{r-fi)U{fi). 

Here the kernel G{r) obeys the equation 

J df2G{f — r2)C{r2 — ^i) = 6{r — fi). (8) 

To perform the actual integration we use the following presentation of the delta function 

6{x) = ^ [ dye^y^, (9) 



27r , 

and then the Gaussian structure of the action S allows us to calculate the integral ([8 

V ( uiv 



K = j dfdfifv{r)C{r- 

This exact result is valid for any Gaussian field U. By definition, C(0) = (f/^) and in 
a typical case C{r) = a^f[r). Hence, K cc and for this reason we introduce a new 
parameter k 

K = Ka\ (11) 

which depends only on the spatial decay of C{r). In future we will omit the factor o"^ in 
C{f). Using the Fourier transforms of C{r) and /y(r*) we obtain 

= ^ 1 dkMk)C{k)M-k), (12) 

and for a spherical domain with radius Ro 

^ Att 

fv{k) = yj (sin kRo — kRo cos kRo) . (13) 

rh 
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A. Noncorrelated field 



Let us analyze Eq. (fTOj) for some particular cases. If U{r) is a field without spatial 
correlations, then 

C{f) = a^6{f), 

and 

K = aV, Pv{Uo) ■■ 



(14) 



V 



1/2 



Note, that in this particular case Py(f/o) depends only on the volume V of the domain and 
not on its shape, as it should be for a totally noncorrelated field distribution. Note also that 
V/ is actually the number s of lattice sites in the domain, so the leading asymptotics for 
s > 1 is 

InP, oc-s. (16) 
The noncorrelated Gaussian random energy is the base of the famous Gaussian disorder 



model, developed 
organic materials 



jy H. Bassler for the description of charge carrier transport in disordered 



SOj. The correlated model 2l|, |22] could be considered as a natural exten- 



sion of the Bassler's model in order to explain the experimental mobility field dependence 
Eq. (ED. 



B. Dipolar-like field in a spherical domain 

Now let us discuss the most interesting case of the dipolar-like correlated field C(r) (x a/r. 
In this case k depends not only on the total volume of the domain, but also on its geometry. 
Using the Fourier transforms of C{r) 

c(M ^ (17) 



we obtain for a spherical domain 

= 2^AaRl (18) 

This result demonstrates a tremendous difference with the noncorrelated case, described by 
Eq. ( |T5|) . because in the leading asymptotics In Py oc —Ro, i.e. it is proportional to the 
linear size of the domain, and not to its volume. The comparison between the analytic 
result (fT8|) and the simulation data is shown in Figs. [2] and [31 Statistics for all figures have 
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FIG. 2: Comparison between simulation data (+,□) and analytic results (llOp and (llSp for the 
dipolar glass with Rq = 5a (+) and Rq = 10a (□), respectively. Note that there are no adjustable 
parameters in this plot. 

been gathered for a basic sample with a size of 256 x 256 x 256 lattice sites with periodic 
boundary conditions and 10,000 realizations of the random field U (apart from Fig. [2], where 
1,000 realizations of the random field were used). Particular distributions of U{r) have been 
generated in the usual way. There is no correlation among the fluctuations in momentum 
space U{k) for different k, so we generated distributions of U{k) and then calculated Fourier 
transform to get U{f). 

If a domain has an arbitrary shape but still could be characterized by a single linear scale 
Rq, then 

KocRl (19) 

just because of dimensionality argument, though the coefficient of proportionality depends 
on the actual shape of the domain. One can rewrite the relation (fT9|) in the following form 

K = gAaV''!^, (20) 

where the coefficient g depends on the shape of the domain and for a sphere = 2(367r)^/^/5. 
The calculation of the coefficient g for a more general case of elliptic domains is presented 
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FIG. 3: Dependence of k on Rq for a dipolar glass. The slope of the straight line equals to 0.95. 
According to Eq. ()18|) . it should be equal to 0.91. The reason for a nonzero intercept is the finite 
size of the basic cell. 

in the Appendix. This calculation shows that g attains a maximum g = go for a spherical 
shape and is significantly smaller than go only for very elongated or oblate ellipsoids. 

III. ESTIMATION FOR CLUSTER NUMBERS 

The number of spherical domains ny(?7o) per unit volume, having an average energy 
greater than Uq, is approximately equal to 



here the coefficient 1 /V reflects the number of non-overlapping independent domains in any 
finite sample. If Uq ^ a, then 



We may expect that Eq. (I22l) gives a reasonable estimation for the number Ug of the true 
clusters, i.e. domains, where U{r) > Uq everywhere (assuming V = a^s), at least for the 
leading term of the asymptotic dependence of Ug on s (the very use of the continuous model 






(22) 
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of the random medium suggests that our consideration is vahd only for s 3> 1). In addi- 
tion, because we consider the distribution of the average field in the most compact domain (a 
sphere), this estimation could be valid only for clusters far away from the percolation thresh- 
old (this is equivalent to Uq ^ a). At the percolation threshold clusters typically have a 



fractal-like structure 



If this assumption is true, then for the noncorrelated Gaussian field 

Ti^ """"P ("^"^^^) ' ^^^^ 

and for the dipolar-like Gaussian field 

oc ^ exp (-5.^^^/^) , (24) 

where we take into account the possibility that for true clusters the coefficients Bnc and 5^ 
might differ from the corresponding values and B^, estimated from Eqs. ( fT5l) and ( fTSj) 
for spherical domains 

SL = I, (25) 

One can reasonably assume that Ba does not differ significantly from because the spher- 
ical domains are the most probable ones (see Appendix). We compared Eqs. (!23|) and (!24|l 
with the simulation data and found that they provide good approximations for the true 
cluster numbers (see Figs. [Hand [5]). 

In order to understand the true status of Eqs. ( I23l) -( l26l) . let us compare the result for the 
noncorrelated field with Eq. ([T]), which is an exact result for the noncorrelated percolation. 
The noncorrelated Gaussian field problem is exactly equivalent to the classic percolation 
problem [l| with 

Comparing Eqs. ([1]), (!23l) . and (p7j) we see that our simple estimation (l23l) provides at least 
the right leading asymptotics for Us 

l„„. = _^,H.„(£|,). (28) 

(note that Bnc = B^c — 1/2); so both the functional kind of the asymptotic dependence of 
Us on s and the coefficient of proportionality are true for the noncorrelated field. In fact, 
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p=;erfc ^ U^-;^exp , f/o/a » 1. (27) 




FIG. 4: Cluster numbers Ug for a noncorrelated Gaussian field. Threshold energy Uo/a varies from 
2.0 to 3.5 (with the step 0.25) from the topmost curve downwards, and the lines are provided as 
guides for an eye (inset). In proper coordinates all curves approximately collapse to a uniform 
straight line with the slope —0.47. According to Eq. p5|) . the slope should be equal to —1/2. 

even the small difference between 1/2 and the corresponding fitting coefficient in Fig. H] 

could be perfectly well explained by the contribution of higher order terms in Eq. ( l28l) for 

f/o/c — 2 — 3. If we fit only the data for sufficiently large values of the threshold energy, 

then the slope becomes even more close to —1/2. 

If we consider the case of a dipolar-like field, than again the fit of the true simulated Ug 

to Eq. (12^ gives Bfi ^ 0.31 which is very close to B'^ ^ 0.34. Again, if we try to fit only 

data points for t/o/cr ~ 3.5 — 4, then the agreement between and 5^ becomes better. We 

would like to put forward the hypothesis that for the DG model the asymptotic expansion 

/rr2 \ 7-/2 

lnn, = -S°^sV3 + J^.Vsj^ ^s'f^^l (29) 

is valid too. If so, we may suggest that the corresponding asymptotics for and s ^ 1, 
f/o ^ cr is valid for any Gaussian field with k, calculated by Eq. f[T^ for a spherical domain. 
This strong hypothesis certainly should be tested more thoroughly, but, nonetheless, our 
simulation data provides important arguments in its favor. Another interesting question is 
how valid are power law corrections to the leading exponents in Eqs. (I23l) and (|2^ . For 
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FIG. 5: Cluster numbers Ug for the dipolar-like Gaussian field. Threshold energy Uo/a varies from 
2.75 to 4.0 (with the step 0.25) from the topmost curve downwards, and the lines are provided 
as guides for an eye (inset). Again, as in Fig. U in proper coordinates all curves approximately 
collapse to the uniform straight line with the slope —0.31. According to Eq. (jlSp . the slope should 
be equal to —0.34. 

the noncorrelated percolation in Eq. ([1]) in the 3D case the exact result is ^ = 3/2 0] and 
agrees with Eq. (l23l) . though this agreement, quite possibly, is an accident. 

If we consider Eq. ( 128|) . it is obvious that it is universal and does not depend on the 
particular structure of the lattice. This is not so, seemingly, for Eq. (l29l) . where the co- 



or different lattices (the 



261]). At the same time. 



efficient depends on parameter A, which, in turn, is different 
particular value A ^ 0.76 is valid only for a simple cubic lattice 
we cannot expect this kind of dependence for s ^ 1, where the particular structure of the 
lattice should be unimportant. This seeming contradiction could be resolved if we recall that 
for the DG model the parameter depends on the lattice too. In fact, the combination 
Aa'^a is invariant 

= ' (30) 



where c is the concentration of dipoles 



261 . Clearly, in such a case the correlation function 



[IH]) does not depend on any microscopic characteristic of the random dipolar medium, while 
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the combination B^/a"^ depends only on the lattice scale a and not on the particular structure 
of the lattice. 

It was found previously that in the correlating percolation problem many features of the 
percolation near the percolation threshold are not, in fact, very sensitive to the correlation. 
For example, in some cases the percolation threshold is the same for correlated and non- 
correlated problems [l^, and cluster numbers sometimes are the same as well 31 1. From 
this point of view it is very interesting that the asymptotic behavior of Ug for deep clusters 
differs significantly for correlated and noncorrelated Gaussian fields. 



IV. CONCLUSION 



In this paper we discussed the distribution of an average value in a finite domain for 
different Gaussian random fields. We found that for very different types of Gaussian fields 
(in terms of their spatial correlation properties) the distribution of the average energy could 
serve as a good estimation for the true cluster numbers per lattice site for large "deep" 
clusters, where s ^ 1 and the threshold energy Uq is significantly greater than the rms 
disorder a. Comparison of the analytical results for Pv{Uo), calculated for a spherical 
domain, and computer simulation data for supports the hypothesis that Pv{Uo) provides 
the exact leading asymptotic term for Ug. In our consideration we discussed particular 
Gaussian fields, relevant to the description of charge carrier transport in disordered organic 
materials. Nonetheless, the suggested approach could be used for other random Gaussian 
fields as well. Generalization to other spatial dimensions (beyond 3D) is also possible. 
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FIG. 6: Coefficient 5e(a) for the ellipsoidal domains; go = 2(367r)^/^/5 is the corresponding coeffi- 
cient for a sphere. Broken lines correspond to approximations (lA.Sh and (|A.6p . 

APPENDIX: DIPOLAR-LIKE FIELD IN AN ELLIPSOIDAL DOMAIN 



Let us consider domains of a non-spherical shape in order to estimate the influence of 
the domain shape on the probability to have a particular value of Uq in such domain. The 
simplest choice is to study the distribution of the fleld in ellipsoidal domains with half-axes 
biRo, 62-R05 and &3-R05 where bi are scale coefficients. A direct calculation shows that for 
ellipsoidal domains 



{2'ny 



dkfvMCe{k)fvo{-k), 



(A.l) 



iirAa 



here the function fvo{k) is exactly the same one as the corresponding function for spherical 
domains and Vq = 47ri?Q/3. For this reason 

1 



-Kohih2hI, 



dx 



(A.2) 



1/2 ■ 
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here kq is the corresponding value for a sphere. Let us calculate this integral for ellipsoidal 
domains having rotational symmetry with bi = 62- In this case 

dx 



^ = bi ^ , ,,,, = (A.3) 



2hl I arctgVa^ _ 1^ « > 1, 



' 2 i_v/TZ^' " ^ J-, 

where a = hi/h^. Taking into account that the volume of the ellipsoidal domain is equal to 
V = 616263 Vo, we obtain 



ge{a) a^/^ J arctgVa^ - 1, a > 1, 

(A.4) 



90 j ij^ii^^ 

and in the limiting cases 



fl'o a 



9e(a) ^ / . 

Eqs. f lA.5l) and flA.6l) mean that domains which differ significantly from the spherical ones 
have much smaller probability to occur (for the same values of Uo and V). The general 
behavior oi g^ia) is shown in Fig. [61 
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